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We present a microscopic model of a charge carrier transfer under an action of a con- 
stant electric field in a complex medium. Generalizing previous theoretical approaches, 
we model the dynamical environment hindering the carrier motion by dynamic perco- 
lation, 1.6. « clS db medium comprising particles which move randomly on a simple cubic 
lattice, constrained by hard-core exclusion, and may spontaneously annihilate and re- 
appear at some prescribed rates. We determine analytically the density profiles of the 
"environment" particles, as seen from the stationary moving charge carrier, and calculate 
its terminal velocity, V c , as the function of the applied field and other system parameters. 
We realize that for sufficiently small external fields the force exerted on the carrier by 
the "environment" particles shows a viscous- like behavior and define an analog of the 
Stokes formula for such dynamic percolative environments. The corresponding friction 
coefficient is also derived. We show that the density profile of the environment particles 
is strongly inhomogeneous: In front of the stationary moving charge carrier the density is 
higher than the average density, p s , and approaches it as an exponential function of the 
distance from the carrier. Behind the carrier the local density is lower than p s and the 
approach towards p s may proceed differently depending on whether the particles number 
is or is not explicitly conserved. 

1. INTRODUCTION 

The percolation concept has been a key notion for understanding transport and con- 
duction processes in a wide range of complex disordered media. A few stray examples 
are ionic conduction in polymeric, amorphous or glassy ceramic electrolytes, diffusion in 
biological tissues and permeability of disordered membranes jll2l3l4j . 

Most of physical situations studied so far (see Refs. |l|2|3|4j and references therein) 
concern systems in which complex disordered environment can be considered as " frozen" ; 
that is, the random environment in which a given transfer process takes place does not 
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change in time. This is certainly the case in many instances, but it is not true in general. 
In fact, there are many experimental systems in which the static percolation picture is not 
valid since the structure of the host material reorganizes itself on a time scale comparable 
to that at which the transfer itself occurs. Such a reorganization happens, namely, in 
certain biomembranes [5 , solid protonic conductors [6], oil-continuous microemulsions 
|7|8|9|lUj or polymer electrolytes |11|12|13] . 

In particular, ionic transport across a biomembrane, such as, e.g., gramicidin- A, pro- 
ceeds by the motion of ions through molecular channels along which they encounter po- 
tential barriers whose heights fluctuate randomly in time. The fluctuations of potential 
barriers may impede the transport constituting an important controlling factor [Sj. In 
the case of protonic conduction by the Grotthus mechanism |BJ, site-to-site hopping of 
charge carriers takes place between neighboring H 2 or NH 3 groups that have a favorable 
relative orientation. Here, the structural host-reorganization process interacting with the 
carrier motion occurs due to thermally activated rotation of the H 2 or NH 3 groups. In 
a similar fashion, within oil-continuous microemulsions, the charge transport proceeds by 
charge being transferred from one water globule to another, as globules approach each 
other in their Brownian motion |7l8l9llUj . Last but not least, in polymer electrolytes, such 
as, e.g. polyethylene oxide complexed non-stoichiometrically with the ionic salt NaSCN, 
the Na + ions are largely tetrahedrally coordinated by polyether oxygens, but at the same 
time that Na + ions hop from one fourfold coordination site to another, the oxygens 
themselves, along with the polymeric backbone, undergo large-amplitude wagging and 
even diffusive motion |ll|12|13j . 

Clearly, all these examples involve two characteristic time scales: the typical time r 
between two successive hops of the charge carrier and a typical renewal time r* of the 
environment itself; that is, the time at which the host complex medium re-organizes itself 
sufficiently enough to provide a new set of available pathways for charge carrier transfer. 
Accordingly, the conventional static percolation picture can be strictly valid when only the 
characteristic time r* gets infinitely large. Only in this case one may expect an anomalous 
large-scale dynamics. On contrary, when r* is finite, dynamic percolation picture has to 
be applied. In this situation, one encounters quite a different behavior when compared 
to the random environments with quenched disorder. As a result, one expects an Ohmic- 
type or Stokes-type linear velocity-force relation for the carrier's terminal velocities as a 
function of the applied field, in contrast to the threshold behavior and anomalous dynamics 
predicted by the conventional static percolation theory. 

The prefactor in the linear velocity-force relation may depend, however, in a non-trivial 
way on the system's parameters and this dependence constitutes the main challenge for 
the theoretical analysis here. On the other hand, we note that in the above mentioned 
examples of the dynamic percolative environments quite different physical processes are 
responsible for the time evolution of the host medium. Consequently, one expects that 
the prefactor in the Stokes-type velocity-force relation should also be dependent on the 
precise mechanism which underlies the temporal re-organization of the environment. 

Theoretical modeling of charge carrier transfer in dynamic percolative environments 
has followed basically two different lines of thought. Early models of dynamic percolation 
|14|15j described the random environment within the framework of a standard bond- 
percolation model allowing the strength of each bond to fluctuate in time between zero 



3 



and some finite value. In this approach, the dynamics of the host medium |14J15j was 
accounted for by a series of instantaneous renewal events. These events were assumed to 
occur at random times, chosen from a renewal time distribution. In the renewal process the 
positions of all unblocked bonds are being reassigned, such that after each renewal event 
a carrier sees a newly defined network. This approach is thus characterized by a global 
dynamical disorder without global conservation laws and correlations, since the entire 
set of random hopping rates is renewed independently of the previous history. Another 
model characterized by a local dynamical disorder has been proposed in Refs. [Tol and 
|17j . and subsequently generalized to the non-Markovian case in Ref. j!8j. The difference 
of this model in regard to the previous one is that here the hopping rates at different 
sites fluctuate independently of each other. That is, states of individual bonds, rather 
than that of the whole lattice change in the renewal events. To describe the dynamical 
behavior in the local dynamical disorder dynamical mean-field theory has been 

proposed |16|17j . based on the effective medium approximation introduced for the analysis 
of random walks on lattices with static disorder [T^]. Subsequently, this model has been 
generalized to include the possibility of multistate transformations of the dynamically 
random complex medium [20]. More recently, several exactly solvable one-dimensional 
models with global and local dynamical disorder have been discussed [2T] . 

The second approach to the problem emerged within the context of ionic conductivity 
in superionic solids. Here, the dynamical percolative environment has been considered as 
a multicomponent mixture of mobile species in which one or several neutral components 
block the carrier component [221 • in particular, such a situation can be observed in a supe- 
rionic conductor /^"-alumina, doped with two different ionic species (e.g. Na + and Ba 2+ ), 
where small Na + ions are rather mobile, while the larger Ba 2+ ions move essentially 
slower and temporarily block the Na + ions. Contrary to the previous line of thought, 
the dynamics of such a percolative environment has essential correlations, generated by 
hard-core exclusion interactions between the species involved. Moreover, a salient feature 
of these situations is that here the total number of the particles involved is conserved, 
i.e., dynamics of the environment is globally constrained by a certain "conservation law". 
The impact of this property on charge carrier transfer will become clearer as we proceed. 

Further on, in Ref. [22] the frequency-dependent ionic conductivity of the light species 
has been analyzed combining a continuous time random walk approach for the dynamical 
problem with an effective medium approximation describing the frozen environment of 
slow species. Next, as an explanation of the sharp increase of electrical conductivity 
transition in water-in-oil microemulsions when the volume fraction of water is increased 
towards a certain threshold value, in Refs. [H] and jU] it has been proposed that the charge 
carriers are not trapped in the finite water clusters, but rather a charge on a water globule 
can propagate by either hopping to a neighboring globule, when they approach each other, 
or via the diffusion of the host globule itself. This picture has been interpreted in terms of 
a model similar to that employed in Ref. (221 > with the only difference being that here the 
" blockers" of Ref. [22] play the role of the transient charge carriers. In the model of Refs. jHj 
and P] , in which the host dynamics is influenced by spatial correlations and conservation 
of the number of the water globules involved, the conductivity depends on the rate of 
cluster rearrangement. Lastly, a similar problem of a carrier diffusion in an environment 
created by mobile hard-core lattice-gas particles has been analyzed in Ref. [23] by using 
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the developed dynamic bond percolation theory of Refs.jEj and [To] . 

In this paper, we pursue our previous analysis of the charge carrier transfer in complex 
media [21] considering a generalize dynamic percolation model which shares common fea- 
tures with both bond-fluctuating models of Refs. |14I15|16|17|18|2U|21] and those involving 
mobile blockers of Refs. [22023] . The system we consider consists of a host lattice, which 
here is a regular cubic lattice whose sites support at most a single occupancy, hard-core 
"environment" particles, and a single, hard-core charge carrier particle subject to external 
electric field. The "environment" particles move on the lattice by performing a random 
hopping between the neighboring lattice sites, which is constrained by the hard-core inter- 
actions, and may disappear from and re-appear (renewal processes) on the empty sites of 
the lattice with some prescribed rates. Contrary to the environment particles the charge 
carrier can not disappear spontaneously, and is subject to a constant external force E. 
Hence, the carrier performs a biased random walk, which is constrained by the hard-core 
interactions with the "environment" particles. In consequence, on may consider the car- 
rier as some "probe" designed to measure the response of the percolative environment to 
the internal perturbation or, in other words, the frictional properties of such a dynamical 
environment. 

Now, a salient feature of our model, which makes it different to the previously pro- 
posed models of dynamic percolation, is that here the interaction between "environment" 
particles and the carrier are included, such that the latter may influence the dynamics 
of the environment. This results, as we proceed to show, in the emergence of compli- 
cated density profiles of the "environment" particles around the carrier. These profiles, 
as well as the terminal velocity V c of the carrier, are determined here explicitly, in terms 
of an approximate approach of Ref . [23] . which is based on the decoupling of the triple 
carrier-particle-particle correlation functions into the product of pair correlations. 

We realize that the "environment" particles tend to accumulate in front of the driven 
charge carrier creating a "traffic jam", which impedes its motion. The density profiles 
around the carrier are asymmetric: the local density of the "environment" particles in 
front of the carrier is higher than the average and approaches the average value as an 
exponential function of the distance from the carrier. The characteristic length and the 
amplitude of the density relaxation function are calculated explicitly. On the other hand, 
behind the carrier, i.e. in its "wake", the local density is lower than the average. We find 
that, surprisingly, the functional form of the density profile might be very different de- 
pending on the condition whether the number of particles in the percolative environment 
is explicitly conserved or not; as a matter of fact, the local density behind the carrier 
may tend to the average value either as an exponential or even as an algebraic function 
of the distance, revealing in the latter case especially strong memory effects and strong 
correlations between the particle distribution in the environment and the carrier position. 
Further on, we find that the terminal velocity of the carrier particle depends explicitly 
on the excess density in the "jammed" region in front of the carrier, as well as on the 
"environment" particles density past the carrier. Both, in turn, are dependent on the 
magnitude of the velocity, as well as on the rate of the renewal processes and the rate 
at which the "environment" particles can diffuse away from the carrier. The interplay 
between the jamming effect of the environment, produced by the carrier particle, and 
the rate of its homogenization due to diffusive smoothening and renewal processes, man- 
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ifests itself as a medium- induced frictional force exerted on the carrier, whose magnitude 
depends on the carrier velocity. As a consequence of such a non-linear coupling, in the 
general case, (i.e. for arbitrary rates of the renewal and diffusive processes), V c can be 
found only implicitly, as the solution of a non-linear equation relating V c to the system 
parameters. This equation simplifies considerably in the limit of small applied external 
fields E and we find that the force-velocity relation to the field becomes linear. This im- 
plies that the frictional force exerted on the carrier particle by the environment is viscous. 
This linear force-velocity relation can be therefore interpreted as the analog of the Stokes 
formula for the dynamic percolative environment under study; in this case, the carrier 
velocity is calculated explicitly as well as the corresponding friction coefficient. Assuming 
that the Einstein relation between the carrier mobility and its diffusion coefficient holds 
we then estimate the self-diffusion coefficient of the carrier in absence of external 
field. We show that when only diffusive re-arrangement of the percolative environment 
is allowed, while the renewal processes are suppressed, the general expression for the dif- 
fusion coefficient reduces to the one obtained previously in Refs.j2Zj and [2H1- We note 
that the result of Refs. [27] and [2E| is known to serve as a very good approximation for 
the self-diffusion coefficient in hard-core lattice-gases [29J. 

The paper is structured as follows: In Section II we formulate the model and introduce 
basic notations. In Section III we write down the dynamical equations which govern the 
time evolution of the "environment" particles and of the carrier. Section IV is devoted 
to the analytical solution of these evolution equations in the limit t — > oo. In Section V 
we derive explicit asymptotic results for the carrier terminal velocity in the limit of small 
applied external fields E and obtain the analog of the Stokes formula for such a percolative 
environment. Asymptotic behavior of the density profiles of the "environment" particles 
around the carrier is discussed in Section VI. Finally, we conclude in Section VII with a 
brief summary and discussion of our results. 

2. THE MODEL. 

The model for charge carrier transfer in complex media consists of a three-dimensional 
simple cubic lattice of spacing a, the sites of which are partially occupied by identical 
hard-core "environment" particles and a single, hard-core, carrier particle. For both types 
of particles the hard-core interactions prevent multiple occupancy of the lattice sites; that 
is, no two "environment" particles or the "carrier" and an "environment" particle can 
occupy simultaneously the same site, and particles can not pass through each other. 

The occupation of the lattice sites by the "environment" particles is characterized by the 
time-dependent occupation variable 77(f), f being the lattice- vector of the site in question. 
This variable assumes two values: 



Next, the dynamics of the "environment" particles is defined via the following rules The 
particles may, at a given rate, spontaneously disappear from the lattice, and may re- 
appear at random positions and random time moments, which is reminiscent of the host 
medium dynamics stipulated in Refs. |14|15|16|17|18l2U|21j . These two processes will be 
referred to, in general, as renewal processes. 




1, if the site f is occupied 
0, if the site f is empty 



(1) 
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Now, the environment particles move randomly by performing nearest-neighbor random 
walks constrained by the hard-core interactions, which is the main feature of the approach 
in Refs. (22I2S1- More specifically, we stipulate that any of the "environment" particles 
waits a time 5t, which has an exponential probability distribution with a mean r*, and 
then chooses from a few possibilities: (a) disappearing from the lattice at rate g/r*, 
which is realized instantaneously, or (b) attempting to hop, at rate l/Qr*, onto one of 6 
neighboring sites. The hop is actually fulfilled if the target site is not occupied at this 
time moment by any other particle; otherwise, the particle attempting to hop remains at 
its initial position, and (c) particles may re-appear on any vacant lattice site with rate 
f/r*. 

Note that, for simplicity, we assumed that the characteristic diffusion time and the 
renewal times of the "environment" particles are equal to each other. These times, i.e. 
T^f, mean creation time r cr and mean annihilation time r an may, however, be different, 
and can be restored in our final results by a mere replacement I — > lr*/r di f, f — > fr* jr cr 
and g -> gr*/r an . 

Note also that the number of particles is not explicitly conserved in such a dynami- 
cal model of the environment due to the renewal processes. However, in the absence of 
attractive particle-particle interactions and external perturbations, the particles distribu- 
tion on the lattice is uniform and the average occupation p(t) = 77(f) of the lattice tends, 
as t — > oo, to a constant value, p s = f /(f + g). This relation can be thought of as the 
Langmuir adsorption isotherm [50] . 

Hence, the limit r^/ — > oo (or, I — > 0) corresponds to the ordinary site percolation 
model with immobile blocked sites. The limit f,g—>0, (r cr , r an — > oo), while keeping the 
ratio f/g fixed, f/g = p s /(l — p s ), corresponds to the usual hard-core lattice-gas with the 
conserved particles number. 

At time t = we introduce at the origin of the lattice the charge carrier, whose role is 
to probe the response of the environment modeled by dynamic percolation to an external 
perturbation. We stipulate that only the carrier out of all participating particles can not 
disappear from the system, and moreover, its motion is biased by some external constant 
force. As a physical realization, we envisage that the carrier is charged, while all other 
particles are neutral, and the system is exposed to constant external electric field E. 

The dynamics of the carrier particle is defined as follows: We suppose that the waiting 
time between successive jumps of the carrier has also an exponential distribution with a 
mean value r, which may in general be different from the corresponding waiting time of 
the environment particles. Attempting to hop, the carrier first chooses a hop direction 
with probabilities 



p M = exp 



where (3 is the reciprocal temperature, e v (or e M ) stand for six unit lattice vectors, u, p = 
{±1, ±2, ±3}, connecting the carrier position with 6 neighboring lattice sites, and (E-e u ) 
denotes the scalar product. We adopt the convention that ±1 corresponds to ±A, ±2 
corresponds to ±Y while ±3 stands for ±Z. The jump is actually fulfilled when the 
target lattice site is vacant. Otherwise, the carrier remains at its position. For simplicity 
we assume in what follows that the external field is oriented along the X-axis in the 
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positive direction, such that E = (E, 0, 0). Note also that for the choice of the transition 
probabilities as in Eq.(J2J), the detailed balance is naturally preserved. 

3. EVOLUTION EQUATIONS. 

We proceed by writing the evolution equations describing the dynamics of the system. 
Let P(R C , rj; t) denote the joint probability that at time moment t the charge carrier 
occupies position R c and all "environment" particles are in configuration rj = {77(f)}. 
Next, let rf^ denote particles' configuration obtained from 77 by exchanging the occupation 
variables of the sites f and r + ep, i.e. 77(f) <-> rj(f + ep), and ff be the configuration 
obtained from 77 by changing the occupation of the site r as 77(f) <— > 1 — 77(f). Clearly, the 
first type of process appears due to random hops of the "environment" particles, while the 
second one stems from the renewal processes, i.e. random creation and annihilation of the 
"environment" particles. Then, summing up all possible events which can result in the 
configuration (R c , 77) or change this configuration for any other, we find that the temporal 
evolution of the system under study is governed by the following master equation: 

d t p(R c , m t) = -^Yl E { p (Rc, »p; t) - p(r c , m t)} + 

+ - E Pa. { - v(Rc)) P{Rc - en, 77; t) - (l - V (R C + e») P(R C , 77; t)} + 

+ ^E {{i-v(^)P{Rc,f;t)-v(^P(R c ,v;t)} + 

+ ^E {v&p(Rc, tf;t)-(i- 77(f)) p(r C} m t)}. (3) 

Note that the terms in the first (resp. second) line of Eq.© describe random hopping 
motion of the "environment" particles (resp. biased motion of the carrier) in terms of the 
Kawasaki-type particle-vacancy exchanges, while the terms in the third and the fourth 
lines account for the Glauber- type decay and creation of the "environment" particles. 

3.1. Mean velocity of the charge carrier and correlation functions. 

The velocity of the charge carrier can be now readily determined from Eq.(j2J). To do 
it, we multiply both sides of Eq.(j2J) by (R c ■ ei) and sum over all possible configurations 
(R c , rj). This yields the following expression for the carrier mean velocity V c (t): 

V c (t) = - {pi(l - k{e x] t)) - p_ x (l - fc(eL i; *))} , (4) 

T 

where k(X;t) denotes the carrier-" environment" particles pair correlation function: 

k(X;t) = J2r ] (R c + X)P(R c , V ;t). (5) 

Rc,v 

Note that k(X;t) is just the density distribution of the "environment" particles, as seen 
by someone residing on the carrier which moves with velocity V c (t). 
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Note now that V c (t) depends explicitly on the local density of the "environment" par- 
ticles in the immediate vicinity of the carrier. Note also that if the "environment" is 
perfectly homogeneous, i.e., k(X; t) = p s , one has that rj(R c + A) and P(R C , ??; t) in Eq.(j5J) 
decouple, which entails, in virtue of from Eq.(@]), a trivial mean-field- type result 



K (0) = (pi-p-i)(i-p s )-- (6) 



The latter (trivial) expression implies that for the perfectly homogeneous, well-stirred 
environment the frequency of jumps of the carrier particles (t -1 ) merely gets renormalized 
by a factor 1 — p s , which gives the fraction of successful jumps. 

The situation, however, appears to be more complex and there emerge essential backflow 
effects: As a matter of fact, the carrier effectively perturbs the spatial distribution of the 
"environment" particles so that stationary density profiles emerge. This can be contrasted 
to the earlier dynamic percolation models |14pi5fl6|17pi8p20|2ip22p28j in which the carrier 
had no impact on the embedding medium and hence there was no re- arrangement of the 
host medium around the carrier particle. As a consequence, in our model k(X;t) ^ p s , 
and k(X; t) approaches p s only at infinite separations from the carrier, i.e. when |A| — > oo. 
Therefore, we rewrite Eq.(J3} in the form 

V c (t) = V™ " -{Pi(*(ei;t) -ps) -P-i(p a - fc(e_i;t))}, (7) 

T 

which shows explicitly the deviation of the mean velocity of the carrier from the mean- 
field- type result in Eq.© due to the formation of the density profiles. 

3.2. Evolution equations of the pair correlation functions. 

Equation (j3J) signifies that in order to obtain V c (t) we have to evaluate k(e±i] t). Multi- 
plying both sides of Eq.flHJ) by T](R C ) and summing over all configurations (R c ,rj), we find 
that k(X; t) obeys 

dtk&t) = ^2(V M -^V_ M )*(X;t)-^i±^*(X;t) + ^ + 

+ -J2T,P^{ 1 -v(Rc + e^))V^(R c + X)P(R c , m t), (8) 
T ^ Rc,v 

where V M denotes the ascending finite difference operator of the form 

V M /(A)=/(A + e M )-/(A), (9) 

and 

^ f 1, if the site r = r > 
r,r 1 0, otherwise. 

The Kroneker-delta term 5r g signifies that the evolution of the pair correlations, Eq.fJSjl . 
proceeds differently at large separations and in the immediate vicinity of the carrier. This 
stems from the asymmetric hopping rules of the carrier particle defined by Eq.(J2J). 

Note next that the contribution in the second line in Eq.fjHJ), which stems out of the 
bias acting on the carrier, is non-linear with respect to the occupation numbers. In 
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consequence, the pair correlation function is effectively coupled to the evolution of the 
third-order correlations. That is, Eq. (jSJ) is not closed with respect to the pair correlations 
but rather represents a first equation in the infinite hierarchy of coupled equations for 
higher-order correlation functions. One faces, therefore, the problem of solving an infinite 
hierarchy of coupled differential equations and needs to resort to an approximate closure 
scheme. 

3.3. Decoupling Approximation. 

Here we employ a simple non-trivial closure approximation, based on the decoupling of 
the third-order correlation functions into the product of pair correlations. More precisely, 
we assume that for A 7^ e v , the third-order correlation fulfil 



The approximate closure in Eq.ffTTj) has been already used for studying related models 
of biased carrier diffusion in hard-core lattice gases and has been shown to provide quite 
an accurate description of both the dynamical and stationary-state behavior. The decou- 
pling in Eq. fjllj) was first introduced in Ref. [22] to determine the properties of a driven 
carrier diffusion in a one-dimensional hard-core lattice gas with a conserved number of 
particles, i.e. without an exchange of particles with the reservoir. Extensive numeri- 
cal simulations performed in Ref.[23] have demonstrated that such a decoupling is quite 
a plausible approximation for the model under study. Moreover, rigorous probabilistic 
analysis of Ref. pTT] has shown that for this model the results based on the decoupling 
scheme in Eq.(|ll|) are exact. Furthermore, the same closure procedure has been recently 
applied to study spreading of a hard-core lattice gas from a reservoir attached to one of 
the lattice sites [32j. Again, a very good agreement between the analytical results and 
the numerical data has been found. Next, the decoupling in Eq. (lll|) has been used in a 
recent analysis of a biased carrier dynamics in a one-dimensional model of an adsorbed 
monolayer in contact with a vapor phase [33], i.e. a one- dimensional version of the model 
to be studied here. Also in this case an excellent agreement has been observed between 
the analytical predictions and the Monte Carlo simulations data 33 j. We now show that 
the approximate closure of the hierarchy of the evolution equations in Eq. fjllj) allows us to 
reproduce in the limit f,g = and f /g = const the results of Refs. |27j and |28j . which are 
known (see e.g. Ref. |29j) to provide a very good approximation for the carrier diffusion 
coefficient in three-dimensional hard-core lattice gases with arbitrary particle density. We 
expect therefore that such a closure scheme will render a plausible description of the car- 
rier dynamics in our three-dimensional generalized dynamic percolation model. We base 
our further analysis on this approximation. 

Making use of Eq. (jll|) . we find from Eq.fjHJ) that the time evolution of pair correlations 
is governed by the following equations. For all A except for A = {0, e±i, e± 2 , e± 3 }, we have 
that k(X; t) obeys 



v(Rc + X)v(Rc + e u )P{R c , n- 1) « k{X; t)k{e v - 1) 




d t k{X; t) 



Lk(X; t) + 




(12) 



6r* 



10 



where the operator L and coefficients A v (t) are given explicitly by 



fx 1 



(13) 



and 



A fl (t) = l + —p tl (l-k{e tl ;t)), 



(14) 



the operator V M being defined previously in Eq.®, /i = {±1,±2,±3}. On the other 
hand, at the sites adjacent to the carrier one has 



where v = {±1,±2,±3}. 

Now, several comments about equations (fl"2|) and (fl~5|) are in order. First of all, let 
us note that Eq. (|15|) represents, from the mathematical point of view, the boundary 
conditions for the general evolution equation (|12)l. imposed on the sites in the immediate 
vicinity of the carrier. Equations (fT2"|) and (fl3j) have a different form since in the immediate 
vicinity of the carrier its asymmetric hopping rules perturb essentially the "environment" 
particles dynamics. Equations ()12)) and (|15|) possess some intrinsic symmetries and hence 
the number of independent parameters can be reduced. Namely, reversing the field, i.e. 
changing E — > —E, leads to the mere replacement of k(ei,t) by fc(e_i;t) but does not 
affect k(e u ; t) with v = {±2, ±3}, which implies that 

k(er,t)(-E) = Jfe(e_i; *)(£), and k(e u ;t)(-E) = k(e u ;t)(E) for z/ = {±2,±3}. (16) 

Besides, since the transition probabilities in Eq.(J2|) obey p 2 = p~2 — P3 — P-3, one 
evidently has that fc(e 2 ;£) = fc(e_ 2 ;t) = k(e 3 ;t) = k(e^ 3 ;t), and hence, by symmetry, 
A 2 (t) = A_ 2 (t) = A 3 (t) = A_ 3 (t), which somewhat simplifies equations (JT^J) and (|T3|). 
Lastly, we note that despite the fact that using the decoupling scheme in Eq.(fTT|) we 
effectively close the system of equations on the level of the pair correlations, the solution 
of Eqs. (|T2^1 and (fTH|) still poses serious technical difficulties. Namely, these equations are 
strongly non-linear with respect to the carrier velocity, which introduces the gradient term 
on the rhs of the evolution equations for the pair correlation, and depends by itself on the 
values of the "environment" particles densities in the immediate vicinity of the carrier. 
Below we discuss a solution to this non-linear problem, focusing on the limit t — > oo. 

4. SOLUTION OF THE DECOUPLED EVOLUTION EQUATIONS IN THE 
STATIONARY STATE. 

Consider the limit t — > oo and suppose that the density profiles and the stationary 
velocity of the carrier have non-trivial stationary values: k(X) = lim^oo k(X; t), V c = 
lim^oo V c (t) and A^ = lim^oo A^{t). As the next step, we define the local deviations of 
k(X) from the unperturbed density: h(X) = k(X) — p s . This yields the following system 
of equations: 




(15) 



Lh{X) = 0, 



(17) 
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for A 7^ {0, e±i, e±2, e±s}, while for the special sites adjacent to the carrier one has 

(L + A u )h(e u )+p s (A u -A^) = 0, (18) 

Equations (|17j) and ()18|) determine the spatial distribution of the deviation from the 
unperturbed density p s in the stationary state. Note also that in virtue of the symmetry 
relations h(e±2) = h(e±s) and A 2 = A_ 2 = A 3 = A_ 3 . 

To solve the coupled non-linear Eqs.(|^|). (fT7|) and (fTH|) we proceed in the following, 
standard manner: We first solve these equations supposing that the carrier stationary 
velocity is a fixed, given parameter (and hence, the functions A u entering Eqs. (jl7j) and (fTSj) 
are known). In doing so, we obtain h(X) in the parameterized form h(X) = h(X; A±x, A 2 ). 
Then, choosing particular values A = {ef±i, e±2, £±3} and making use of the definition of 
Ap, we find a system of three linear equations with three unknowns of the form 

6t* 

A v = l + —pJl-p a -h(e„;A ±u A 2 j), (19) 



h 

where v = {±1,2}, which will allow us to obtain a closure relation and hence, to define 
all A v explicitly (and hence, all h{e v )). Finally, substituting the results into Eq.(HJ), which 
can be written down in terms of A v as 

V C =^(A 1 -A. 1 ), (20) 
we arrive at a closed-form equation determining implicitly the stationary velocity. 

4.1. Density profiles in the dynamic percolative environment. 

The general solution of Eqs. (jl7j) and (fTHj) can be obtained in a standard fashion by 
introducing the following generating function: 

H(w h w 2 ,w 3 )= K\)w^w¥wfi (21) 

where ni,ri2 and n 3 are the components of the vector A, A = ein 1 +e2n 2 +e3n 3 . Multiplying 
both sides of Eqs.()17|) and ()18j) by w^w^w^ 3 and performing summation, we find then 
that H(wi,W2,w 3 ) is given explicitly by 



( £* (A v (w^ v] - l)h(e u ) + Ps (A» - A- U )w 



H(w h w 2 , w 3 ) = -i — v - ' : ' ' " v " (22) 

an expression which allows us to determine the stationary density profiles as seen from 
the carrier which moves with a constant velocity V c . 

Inverting next the generating function, Eq. (}2*2*|) . we get, after rather lengthy but straight- 
forward calculations, the following explicit result for the local deviation from the unper- 
turbed density: 

h(X) = "£A v h(S v )V- v - p s (A 1 - A_0(Vi - V_i)}F(A), (23) 

V 

where F(X) is given by 

*® - &f r <->- (^) >•« (4*) (4*) * <*> 
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and 

A | g(/ + g) , , , i 6(/ + ff) 

a = 2^ A/ H ; = + A_i + 4A 2 H (25) 

„ t t 

Consequently, the particles density distribution as seen from the carrier moving with a 
constant velocity V c obeys 

k(X) = Ps + a-'{ J2 A v h{e v V- v - Ps (A 1 - A_ a )(Vi - V_i)}F(A), (26) 

where we have to determine three yet unknown parameters A x , A-\ and A 2 . 

To determine these parameters, we set in Eq.(J23j) X — e x , X — e_i and A = e 2 , which 
results in the system of three closed-form equations determining the unknown functions 
A„,i/ = {±1,2}, 

A, = 1 + (l - Ps - Ps (A 1 - , (27) 

It { detC J 

where C is a square matrix of the third order defined as 

'A x V- X F(e x )-a A_!ViF(ei) A 2 V_ 2 F(e*i) \ 

A 1 V_ 1 F(e_ 1 ) A_!ViF(e_i) - a A 2 V_ 2j P(e_ x ) (28) 
AiV-xFCe'a) A^ViF^) A 2 V_ 2 F(e 2 ) - a/ 

while stands for the matrix obtained from C by replacing the u-th column by a column 
vector ((Vi — V_i)F(e 1/ )) !/ . Equation (J2EJ), together with the definition of the coefficients 
A v , constitutes the first general result of our analysis defining the density distribution in 
the percolative environment under study. 

4.2. General force-velocity relation. 

Substituting Eq.(|2*7j) into (|2"U|). we find that the stationary velocity of the carrier particle 
is defined implicitly as the solution of equation: 

t/ °( \fi \!i i 6r*pidetCi -p_i detC_i\ 1 , . 

V c =-( Pl -P-i)(l-p.){l + P.- F ^ , (29) 



where C x and C_i are the following square matrices of the third order: 

C-i 



( (Vi - V_i)F(ei) A_ 1 ViF(g , i) A 2 V_ 2 F(e!) 
(Vi - V-i)F(eLi) A. 1 ViF(e_i) - a A 2 V_ 2 F(e_!) I (30) 



V (Vi - V_i)F(e 2 ) A_iViF(e 2 ) A 2 V_ 2 F(e 2 ) - a, 



and 



(A X V- X F{? X ) -a (Vi - V_i)F(ei) AV^F^) 

AiV-xF^O (Vi - V_i)F(eU) A 2 V_ 2 F(e_!) j . (31) 

V AV_!F(e 2 ) (Vi - V_!)F(e 2 ) A 2 V_ 2 F(e 2 ) - a 



Equation (J2"H|) represents our second principal result defining the force-velocity relation 
in the dynamic percolative environment for an arbitrary field and arbitrary rates of the 
diffusive and renewal processes. 
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5. CARRIER VELOCITY, FRICTION AND DIFFUSION COEFFICIENTS. 

Consider now the case when the applied external field E is small. Expanding the 
transition probabilities p, and p_i in the Taylor series up to the first order in powers of 
the external field, i.e., setting 

we find that V c , Eq.flUJ), follows 

V c ~ £{<tPE(1 - p s ) - (h^) - Me-i))}- (33) 
On the other hand, Eq. (j2*3*j) entails that 

where 

A = limA u = 1 + ^(1- p 8 ), (35) 
E->o It 



a = hm a = 6 1 1 H : 1 : — ), (36) 



and 
while 

C(x) = |^°° e~*lg(:rt) (I (xt) - h(xt)) dij 1 = {P(0; 3z) - P(2ei; 3x)} _1 , (37) 
with P(r; £) being the generating function, 

+oo 

P(f;0 = E P iW> (38) 
i=o 

of the probability that a walker starting at the origin and performing a Polya random 
walk on the sites of a three-dimensional cubic lattice will arrive on the j-th step to the 
site with the lattice vector f j3] . 

Consequently, we find that in the limit of a small applied field E the force-velocity 
relation in Eq. (J29|) attains the physically meaningful form of the Stokes formula E = (V c , 
which signifies that the frictional force exerted on the carrier by the environment particles 
is viscous. The effective friction coefficient ( is the sum of two terms, 

C = Co + (coop (39) 

where the first term represents a mean-field-type result (o = 6r//5a 2 (l — p s ) (see Eq.(7)), 
while the second one, (coop, obeys 

c = l lMl (40) 

CO ° P f3aH(l - Ps ).(a C(2A /a ) - A ) 
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The second contribution has a more complicated origin and is associated with the coop- 
erative behavior - formation of a inhomogeneous stationary particle distribution around 
the carrier moving with constant velocity V c . Needless to say, such an effect can not be 
observed within the framework of previous models of dynamic percolation, since there the 
carrier does not influence the host medium dynamics |14I15|16|17|18I20I21I22 23J. 

Let us now compare the relative importance of two contributions, i.e. Co an d (coop, to the 
overall friction. Straightforward analysis shows that the cooperative behavior dominates 
at small and moderate / (which entails also small values of g), while for larger /, when (/(o 
tends to 1, the mean- field behavior becomes most important. The cooperative behavior 
also appears to be more pronounced at larger densities p s . 

Consider next some analytical estimates. We start with the situation, in which diffusion 
of the environment particles is suppressed, i.e. when / = 0. In this case, we get 

(coop _ 2P.S / A I \ 

rj v )-iV (l) 

y 



Co {I - Ps )(lC{y) - I 
where 



+ (42) 



3 V r*(l-p s ) 

Suppose first that p s is small, p s <C 1. Then, y « 1/3(1 + r/r*(f + gj) and we can dis- 
tinguish between two situations: when r <C {f + g)/i~*, i.e. when the carrier moves faster 
than the environment re-organizes itself, and the opposite limit, r ^> (f + g)/T*, when the 
environment changes very rapidly compared to the motion of the carrier. In the former 
case we find that y « 1/3, which yields CWCo « 2p fl /(6£(l/3) - 1), £(1/3) « 0.7942, 
while in the latter case we have y ~ t* /3r(f+g) and ( coop /(o ~ p s T* /3r(/ +g). Note, that 
in both cases the ratio ( CO op/ Co appears to be small, which signifies that at small densities 
p s the mean-field friction dominates. Such a result is not counterintuitive, of course, since 
in the absence of the particles' diffusion, which couples effectively the density evolution 
at different lattice sites, no significant cooperative behavior can emerge at small densities. 
On the other hand, at relatively high densities p s ~ 1 and r/(l — p s ) 3> r* j (/ + g) r, 
when the carrier moves at much faster rate than the host medium reorganizes itself, we 
find that ( CO o P /(o ~ t */3t(/ + g) ^> 1. This result stems from the circumstance that in 
sufficiently dense environments modeled by dynamic percolation a highly inhomogeneous 
density profile emerges even in the absence of particles diffusion. Here, on the one hand, 
the carrier perturbs significantly the particle density in its immediate vicinity. On the 
other hand, the density perturbation created by the carrier does not shift the global bal- 
ance between creation and annihilation events, i.e. the mean particle density still equals 
p s . The latter constraint induces then appearance of essential correlations in particles 
distribution and hence, appearance of cooperative behavior. 

Let us consider the opposite case when the renewal processes are not allowed, which 
means that the particles number is conserved and local density in the percolative envi- 
ronment evolves only due to particles diffusion. In this case we find 



(coop 2t p s 



Co " (/r + r*(l-p s ))(6£(l/3)-l) 



(43) 
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Here, the ratio Ccoop/Co can be large and the "cooperative" friction dominates the mean- 
field one when It <C T*(3p s — 1), which happens, namely, at sufficiently high densities and 
in the limit when the carrier moves at a much faster rate than the environment reorganizes 
itself. Otherwise, the mean-field friction prevails. 

To estimate the carrier particle diffusion coefficient D c we assume the validity of the 
Einstein relation, i.e. j3D c = (see, e.g., Ref.j2E])- We find that, in the general case, 
the carrier diffusion coefficient D c reads 

Dc = ^ j t _ ^ {aoC[2Ao/aa) _ 1 + ^-j (44) 

In the particular case of conserved particles number, when /, g — > but their ratio f/g is 
kept fixed, f/g = p s /(l — p s ), the latter equation reduces to the classical result 

, NK _a 2 (l-p s ) f 1 2p s T* (aA 1 , T *(3p s -1)^ 



D N C K = v 6r |1 - ^-(6A £(l/3) - 1 + ^ ^ j j, (45) 

obtained earlier in Refs.j2TJ and [2E] by different analytical techniques. The result in 
Eq. (j43j) is known to be exact in the limits /) s < 1 and p s ~ 1, and serves as a very 
good approximation for the self-diffusion coefficient in hard-core lattice gases of arbitrary 
density [21]. 

Finally, in the absence of particle diffusion (fluctuating-site percolation), our result for 
the carrier particle diffusion coefficient reduces to 

DT = ^ I 1 - 2P,(4[(1 - Ps) + (/ + 9)r/T*]C(y) + 3p s - lY') (46) 



6r 

Note, however, that this result only applies when both / and g are larger than zero, such 
that the renewal processes take place. In fact, the underlying decoupling scheme is only 
plausible in this case. Similarly to the approximate theories in Refs. [27] and [28], our 
approach predicts that in the absence of the renewal processes D^ er vanishes only when 
p s — > 1, which is an incorrect behavior. 

6. ASYMPTOTIC BEHAVIOR OF THE DENSITY PROFILES. 

The density profiles at large separations in front of and past the carrier can be readily 
deduced from the asymptotical behavior of the following generating function 

+oo 

N(wi)= £ h( ni ,n 2 = 0,n 3 = 0)w^. (47) 

Inversion of Eq. (|22|) with respect to the symmetric coordinates n 2 and n 3 yields then 

, (A/i(ei) + Ps {A 1 - U - l) + (A^hie-i) ~ Ps(A 1 - A_0) (V 1 - l) 

N(w x ) = ± ^ ^ >- 

a — Aiw 1 — A-iWi 

r i T 2^ 2A 2 u , 4A 2 /i(e 2 ) 

x / exp-a;I — ^ x)dx + — - x 

Jo a — A\Wi — A_iu>i a — Aiw 1 — A-\W\ 

f°° 2A 2 ( 2A 2 

/ exp-a;I — - x z)Ui( — =i x)- 

Jo a — AiW 1 — A-iWx v a — AiW 1 — A_iW\ 

2 

a — A\Wi 1 — A-\Wi 



x 



9 A 

- Io( z -i 2 z x))dx (4 



16 



We notice now that N(wi) is a holomorphic function in the region Wi < w\ < W 2 , where 
a - AA 2 ha - AA 2 \2 



Wi 
and 



2A_ X 



2A_ X 



(49) 



_ a-AA 2 //a-4A 2 \2 

Vv2 — — —, \ r 



(50) 



As a consequence, the asymptotic behavior of h(ni,n 2 = 0,ns = 0) in the limit n\ — > 00 
(resp. rii — > —00) is controlled by the behavior of N(wi) in the vicinity of wi = W 2 
(resp. w\ = Wi) (see, for example, the analysis of the generating function singularities 
developed in Ref.|Mj). 

6.1. Asymptotics of the density profiles at large separations in front of the 
carrier. 

Consider first the asymptotic behavior of the density distribution of the "environment" 
particles at large separations in front of the carrier. We find then that in the limit 
Wi — > W 2 , the function N(w-\) follows 



N( Wl ) 



jA 1 h{e 1 )+p s {A 1 -A_ 1 ))[W 2 -l 



Wl— >W2 



A-kA 2 



+ 



4nAo 



In (W 2 - 



(51) 



Then, (cf, Ref. |2]), we obtain the following asymptotical result 



/i(ni,0,0) ~ ni . 



K + 
1 

m 



-ni/A-| 



(52) 



where the characteristic length A + is given explicitly by 



1 ( a/2 - 2A 2 



'a/2 - 2A 2 



(53) 



and the amplitude K + obeys 
k+ _ >i/i(e 1 )+p s (A 1 -A_ 1 ))(W 2 -l 

+ 



+ 



47rA 2 

(A^h(e^) - Ps (Ai ~ A-i) (W, 1 - l) 



>0, 



(54) 



which signifies that the density of the "environment" particles in front of the carrier is 
higher than the average value p s and approaches p s at large separations from the carrier 
as an exponential function of the distance. 



17 



6.2. Asymptotics of the density profiles at large separations behind the carrier. 

We consider next the asymptotic behavior of the "environment" particles density pro- 
files past the carrier particle, which turns out to be very different depending on whether 
the dynamics of the percolative environment obeys the strict conservation of the "envi- 
ronment" particles number or not (the renewal processes are suppressed or allowed). 

6.2.1. Non-conserved particles number. 

In the case when particles may disappear and re-appear on the lattice, one has that 
the root Wi < 1. We find then, following essentially the same lines as in the previous 
subsection, that 



N( Wl ) 



wi— >Wi 



{A l h{ei)+Ps{A l -A_ l )){w l -l 



+ 



AixA 2 

(A^h(e^) - Ps (Ai - ^-i)) (Wf 1 - l) 



+ 



A-kAo 



ln( 



wi - Wi 



(55) 



Hence, in the non-conserved case the approach to the unperturbed value p s is also expo- 
nential when ni — > — oo, and follows 



"711,0,0 ~n r 

where 
A_ = In" 1 



K~ 

Kl 



a/2- 2A 2 



V 



A. 



\ 



'a/2-2A 2 



A 

A., 



(56) 



(57) 



and 
K~ 



{A 1 h{e 1 )+ Ps {A l -A^))(w 1 -l 



A-ih{e. 



4nA 2 
p s {A 1 -A_ 1 



Wf 1 



4ttAo 



< 



(5f 



which implies that the particles density past the carrier is lower than the average. Note 
that, in the general case, A + < A_, which means that the depleted region past the carrier 
is more extended in space than the traffic-jam-like region in front of the carrier. The 
density profiles are therefore asymmetric with respect to the origin, n\ = 0. 

6.2.2. Conserved particles number. 

Finally, we turn to the analysis of the shape of the density profiles of the percolative 
environment behind the carrier in the particular limit when the host medium evolves only 
due to diffusion, while creation and annihilation of particles are completely suppressed. In 
this case, in which the particles number is explicitly conserved, one has that for arbitrary 
value of the field and particles' average density, the root Wi = 1 and, consequently, the 
form of the generating function is qualitatively different from that in Eqs. (j51)) and (|55|) . 



N( Wl ) 



>i+ 



{A x h&) - A_i/i(eLi) 2 Ps {A 1 - A_ x )) 



AttA 9 



(w x - 1) In ( 



(59) 
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Equation (JoUj) implies that in the limit when the particle number is conserved the large-ni 
asymptotic behavior of /t ni) o,o is described by an algebraic function of n\ with a logarithmic 
correction; that is, 

where is an ni-independent constant. Remarkably, the power-law decay of correla- 
tions implies existence of a quasi-long-range order in the percolative environment past the 
carrier. In the conserved case the mixing of the three-dimensional percolative environ- 
ment is not very efficient and there are considerable memory effects - the host medium 
remembers the passage of the carrier on large space and time scales. 

7. CONCLUSIONS 

To conclude, we have presented a microscopic model describing the dynamics of a charge 
carrier, driven by an external field E in a three-dimensional complex medium modeled by 
dynamic percolation, i. e. represented as a cubic lattice partially filled with mobile, hard- 
core "environment" particles which can spontaneously disappear and reappear (renewal 
processes) in the system with some prescribed rates. Our analytical description of the 
transfer process in such a medium has been based on the master equation, describing 
the time evolution of the system, which has allowed us to evaluate a system of coupled 
dynamical equations for the charge carrier velocity and a hierarchy of correlation functions. 
To solve these coupled equations, we have invoked an approximate closure scheme based 
on the decomposition of the third-order correlation functions into a product of pairwise 
correlations, which has been first introduced in Ref.j2S] for a related model of a driven 
carrier dynamics in a one-dimensional lattice gas with conserved particles number. Within 
the framework of this approximation, we have derived a system of coupled, discrete-space 
equations describing evolution of the density profiles of the environment, as seen from the 
moving charge carrier, and its velocity V c . We have shown that V c depends on the density 
of the "environment" particles in front of and past the carrier. Both densities depend 
on the magnitude of the velocity, as well as on the rate of the renewal and diffusive 
processes. As a consequence of such a non-linear coupling, in the general case, (i.e. for 
an arbitrary driving field and arbitrary rates of renewal and diffusive processes), V c has 
been found only implicitly, as the solution of a non-linear equation relating its value to 
the system parameters. This equation, which defines the force-velocity relation for the 
dynamic percolation under study, simplifies considerably in the limit of small applied 
field E. We find that in this limit it attains the physically meaningful form of the Stokes 
formula, which implies, in particular, that the frictional force exerted on the carrier by the 
environment modeled by dynamic percolation is viscous. In this limit, the carrier velocity 
and the friction coefficient are calculated explicitly. In addition, we determined the self- 
diffusion coefficient of the carrier in the absence of the field and show that it reduces to the 
well- know result of Refs. |2Z] and [2H] m the limit when the particles number is conserved. 
Further more, we have found that the density profile around the carrier becomes strongly 
inhomogeneous: the local density of the "environment" particles in front of the carrier is 
higher than the average and approaches the average value as an exponential function of 
the distance from the carrier. On the other hand, behind the carrier the local density is 
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lower than the average, and depending on whether the number of particles is explicitly 
conserved or not, the local density past the carrier may tend to the average value either 
as an exponential or even as an algebraic function of the distance. The latter reveals 
especially strong memory effects and strong correlations between the particle distribution 
in the environment and the carrier position. 
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